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Abstract 

We study thermal rectification in single-walled carbon nanohorns (SWNHs) by using non- 
equilibrium molecular dynamics (MD) method. It is found that the horns with the bigger top 
angles show larger asymmetric heat transport due to the larger structural gradient distribution. 
This kind of gradient behavior can be further adjusted by applying external strain on the SWNHs. 
After being carefully elongated along the axial direction, the thermal rectification in the elongated 
SWNHs can become more obvious than that in undeformed ones. The maximum rectification 
efficiency of SWNHs is much bigger than that of carbon nanotube intramolecular junctions. 

PACS numbers: 66.70. +f, 44.10.+i, 61.46.Fg, 65.80.+n 



electronic mail |wugaxp@gmail.com 
'electronic mail: |phylibw@nus.edu^ 



1 



Although human being has mastered the techniques to control the electronic flux by using 
elements such as diodes and/or transistors for a very long time, how to fabricate such devices 
to control the heat flux is still an open and challenging question. Given the importance of 
heat to human being, thermal controlling devices should also play an important role in our 
future life, ranging from nanoscale calorimeters, cooling of microelectronic processors, to 
macroscopic refrigerators and energy-saving buildings. 

The extensive study of heat conduction in low dimensional systems from microscopic point 
of view in past two decades |l, O, [s] has enriched our understanding of the heat conduction 
in flnite size systems. In turn, some interesting inventions have been made recently. For 
example, two interesting solid state thermal diode/rectifler models have been proposed from 

h n n 

nonlinear lattices. [4, |5|, la] Based on the negative differential thermal resistance (NDTR), 
a thermal transistor model has been constructed. [7] By combining thermal transistors in 
different ways, elementary thermal logic gates such as NOT, AND/OR gates have been also 
demonstrated recently. More significantly, the two segment model of thermal rectifier pro- 
posed in Ref y, has been experimentally realized by using gradual mass-loaded carbon and 
boron nitride nanotubes. It has been illustrated that the carbon nanotube intramolecular 
junctions also exhibits thermal rectifying effect. However, the maximum rectification ef- 
ficiency of the real materials seems to be no more than 1.5. Therefore, to find more efficient 
thermal rectifier is becomingof primary importance for future application. 

Nano carbon materials, HI] such as carbon nanohorn aggregates, 12] have recently at- 
tracted considerable attention since they are promising nanoscale materials with potential 
applications in broad fields, such as catalyst supports in fuel cell electrodes, [3] in gas storage 
devices, Q, Q as well as in drug delivery systems. [3] 

The nanohorn is composed of a horn-shaped, cap-closed carbon nanotube (CNT) formed 
out of a single- walled graphene sheet and several of them are put together to form a fiowerlike 
aggregate with a rather uniform diameter of about 100 nm. CNHs have very large surface 
areas which enable them to support very fine catalyst particles and to entrap drug molecules 
and even gases. 

If the CNH is made of one single graphene sheet, we call it single- walled carbon nanohorn 
(SWNH). Two examples are shown in Fig. 1. The SWNH has an angle of about 19 ° , a 
diameter of 2-4 nm, and a length of about 50 nm. |l^ They come together to form a spherical 
aggregate, and at their center the SWNHs coalesce (see Fig. 6 of Ref. Illl ). 
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In this paper, we will demonstrate that SWNH can be a very efficient thermal rectifier. 
Moreover, we find that the force constants along the SWNHs' axes can be adjusted in 
the elongated structure. With suitable tensile stress, the rectification in deformed SWNH 
becomes more obvious than that in undeformed one, which makes the deformed SWNH an 
ideal candidate for thermal rectifier. 

The SWNHs are obtained by cutting out 60 ° and 120 ° sectors from a graphene sheet 
and then gluing together the two cut sides of the sheet. They are illustrated in Fig. 1 and 
named as SWNH-60 and SWNH-120 in the following text. In addition, we also perform 
simulation on a smaller SWNH-120 structure, the central axial length of which is only half 
of the original SWNH-120. This smaller structure is named SWNH-120-SEG. In Fig. 1, the 
blue atoms are fixed in MD process, and orange atoms are put in thermostats, which are 
realized by the Nose-Hoover thermostat, [l?! The temperatures of left and right ends are Tl 
and Tr, respectively. When ^ Tr, temperature gradient can be built up in the SWNH. 
For convenience, we introduce two variables: AT = ^-|^^, and (T) = '^^^^ . In this work, 
(T) is always kept at 290 K. The heat fiux running from the narrow end to the wide one is 
defined as positive. 

The velocity Verlet method is adopted to integrate the equations of motions with the 
time step of 0.51 fs. The typical total MD process is 1 x 10^ steps which is about 5.1 ns, 
and the statistic averages of interesting quantities start from half of the MD process, i.e., 
5 X 10^ steps are used to relax the system to stationary state. 

The instant temperature of atom i is defined as Tj (t) = ^ [y^ (t)^ + Vy iff + iff), 
where is the mass, v (t) is the velocity at time t and ks the Boltzmann constant. This is 
a result of energy equipartition theorem. 

The thermal fiux is obtained from the thermostats, i.e., the total work from the ther- 
mostats can be regarded as the heat fiux runs from thermostats to the system. For more 
details, please refer to our recent work.[10|] The second-generation reactive empirical bond 
order potential [l8| is employed to describe the C-C interactions. 

We use quantity L/Lq to describe the strength of tensile strain, where L and Lq are 
the central axial length after and before elongation, respectively. Lq is obtained after fully 
structural optimization. 

The rectification, which was defined as heat fiux ratio |J„/Jp| = — Jn/</p,j4, 5] is used to 
quantify the asymmetric heat conduction. Where J„ means the heat fiux runs from wide 
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end to the narrow one, and Jp means the heat flux runs in the opposite direction. 

We first study the heat conduction of SWNHs with different top angles under different 
temperature difference (AT =-0.05, 0.05, -0.30 and 0.30) and/or tensile strains. The results 
are shown in Fig. 2. 

It is obvious that | J„| > \ Jp\ in most cases. That's to say, the heat flux runs preferentially 
along the direction of radius and mass decreasing. This result bears an analogy to the recent 
experimental one, 9| which implies that radius and mass gradients are possible reasons of 
thermal rectification. 

In order to understand the origin of the thermal rectification, we can consider the similar 
phenomena in carbon nanotube intramolecular junctions (IMJs) [10] as a reference. Both of 
the SWNHs and IMJs can be considered as curved graphene surfaces, so their phonon spectra 
are similar to that of the graphene, with some difference induced by the curvature effect and 
Born-von Karman boundary condition along the circumferential direction. For the narrow 
end of the SWNHs, the curvature effect is much significant and the phonon spectra deviates 
much from that of the graphene; whereas for the wide end, the phonon spectra approaches 
to that of the graphene. Therefore, the phonon spectra along the central axis of the SWNH 
changes gradually from narrow to wide end. As we have known, the matching and/or 
mismatching of the energy spectra along direction of the heat transport is the underlying 
mechanism of the rectification in nonlinear lattice systems [4, [s], O] and IMJs,[]I)i] we can 
easily understand the appearance of the thermal rectification in SWNHs. Furthermore, it 
was shown that heat can more easily transport from wider side to narrower one in IMJs. This 
behavior is also consistent with our observation of the thermal rectification in SWNHs. This 
correspondence also implies the common underlying mechanism of the thermal rectification 
in both systems. 

The SWNH-120-SEG can be elongated more than both SWNH-60 and SWNH-120. In 
Fig. 2, the maximum possible tensile strain for SWNH-120-SEG is about 1.06, while that 
for SWNH-60 and SWNH-120 is about 1.03. This is because the largest strain occurs at the 
narrow end, when system is elongated, the narrower end of longer SWNH will endure more 
strain than that of shorter ones due to the higher structural asymmetry. 

It is found from the inset of Fig. 2 that both the temperature difference (|AT|) and the 
tensile strain {L/Lq) have large influence on the heat flux J, but their effects are in contrast 
with each other. Namely, an increase in temperature gradient will cause an increase of heat 



4 



flux, but when tensile strain becomes stronger, heat flux J decreases. When the system is 
elongated, the interatomic force constants decreases, which leads to the decrease of the heat 
flux. 

In addition, more interesting are the influences of the temperature difference and the 
tensile strain on the thermal rectification. 

We find that an increase in temperature difference |AT| will result in an increase in rec- 
tification. The rectification is found to be mainly controlled by the middle-frequency optical 
phonon modes, which can be further excited when temperature gradient increases. In 
fact, these middle-frequency optical phonon modes are most sensitive for curvature effect 
in CNT,[19)] so their participation in transport process also play most important role in the 
thermal rectification of SWNHs. In other word, these phonon modes of narrower end is sen- 
sitive to temperature change, which makes thermal rectification changes when temperature 
gradient changes. 

Under the same |AT|, SWNH-60 shows smallest rectification when the strain is weak. 
This is reasonable because the top angle 9 of the horn can be regarded as an indicator of the 
structural asymmetry and curvature gradient along central axis. When 6' = is considered 
as an extreme situation, the SWNH will degenerate to SWNT and phonon spectra gradient 
disappears, and as we known, there does not exist thermal rectification in SWNT. Moreover, 
the rectification in SWNH-120 is even larger than that in SWNII-120-SEG. This is due to 
the larger structural difference between the two ends of SWNH-120. Therefore, if we want 
to design a device exhibiting large rectification, the best choice is to find a system with 
large structural asymmetry. On the other hand, considering the mass density along the Z 
direction is not homogenous, the SWNHs can be also thought as a mass graded system. 



Recent study 



20| has demonstrated that the thermal rectification exists in mass graded 



anharmonic system. So the mass gradient system, which is a special case of structural 
asymmetry system, is hopefully to be a thermal rectifier. In fact, recent experiment of 
thermal rectifier ^1 has proven this idea. 

As shown in Fig. 2, when tensile strain becomes stronger, the rectification does not 
change monotonically. However, the trend of the rectification versus strain seems similar 
for all the horns. Under weak strain, the rectification increases with the increasing strain; 
nevertheless, the rectification will decrease if the strain becomes large. We can understand 
this as following qualitatively. 
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When the system is elongated, the interatomic force constants do not decrease in the same 
rate along the Z direction. Actually, the end with large radius can hardly be elongated. As 
a result, high local strain occurs near the narrow end and the force constants along the 
Z direction will change gradually. Consequently, the nanohorn becomes a force constant 
graded system. Therefore, the asymmetry in curvature and mass distribution needs to 
compete (or cooperate) with the asymmetry in force constant distribution. Furthermore, 
the force constants are not linearly dependence on the bond length when the strain is strong, 
so the gradient in force constants does not always increase with the increase of strain. Above 
multifold reasons finally lead to complicated change of the phonon spectra gradient, and in 
turn to the complicated trend of the rectification. Some additional details of the trend of 
the curves should due to the different responses of the horns to tensile stress. 

So far, one finds that the deformed SWNH-120-SEG exhibits obvious thermal rectifica- 
tion, and its dependence on strain seems similar to the other two nanohorns in the paper. 
Thus we can take SWNH-120-SEG as a typical example in the following text. In order to find 
the maximum performance of SWNH-120-SEG, we adjust only the temperature difference 
of two ends by fixing LjL^ to 1.015. The results are presented in Fig. 3. For comparison, 
the temperature difference dependence of undeformed SWNH-120-SEG is also shown. 

In the inset of Fig. 3, the heat flux appears obviously asymmetric in both structures, but 
the elongated structure shows larger asymmetry. It can be found that the dependence of heat 
flux on temperature difference rapidly deviate from linear relationship when |AT| is larger 
than about 0.1. And when |Ar| is larger than 0.5, the rectification can be as large as 2 in 
deformed SWNH-120-SEG. The overall maximum rectification is about 2.4 when | AT| ^ 0.8. 
This behavior seems very similar to the junction situation, [lOi] but the rectification in SWNH 
is more obvious. This result is due to more distinct structural asymmetry in SWNHs. 

More significative thing is that a larger temperature difference can induce a smaller heat 
flux when |Ar| > 0.8. This phenomenon is called negative differential thermal resistance 
(NDTR)|3], which is essential for thermal transistor. This indicates that the SWNH might 
be used to construct thermal transistors and other thermal devices such as thermal logic 
gates. 

In summary, thermal rectification, i.e., asymmetric heat conduction, has been studied 
in SWNHs by MD method. The structural asymmetry is shown to have great effect on 
the asymmetric heat transport. The horns with bigger top angle can show larger thermal 
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rectification. Our results also show that SWNHs behave like mass graded systems. 

Furthermore, force constants can be adjusted by external stress to infiuence the structure 
ununiformity. The competition between the radius gradient and force constant gradient 
makes the system show complicated response to the tensile strain. After choosing suitable 
tensile strain, the thermal rectification in the SWNHs may become more obvious than that 
in undeformed ones. 

Our work suggests a possible efficient heat control device, and the efficiency of this device 
can be controlled by simply applying external stress. In this system, on-site potential is not 
requested to achieve asymmetric heat conduction. Because SWNHs has been fabricated 
experimentally, we believe the thermal controlling devices can be made in the near future. 

This work is supported by an Academic Research Grant, R- 144-000-203- 112, from the 
Ministry of Education of the Republic of Singapore, and the DSTA under Project Agreement 
No. POD0410553. 
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FIG. 1: (Color online). The structures of the two SWNHs in our work. Upper panel shows the 
side view and top view of SWNH-60 structure. Lower panel shows the side view and top view of 
SWNH-120 structure. 



Figures 
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FIG. 2: (Color online). Thermal rectification in SWNHs under different tensile strain versus top 
angle and temperature difference. The negative heat fluxes correspond to negative temperature 
differences AT. Inset shows the corresponding heat flux. Error bars are also plotted. 
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FIG. 3: (Color online). Thermal rectification of SWNH-120-SEG under different temperature 
difference. Inset shows the corresponding heat flux. Error bars are also plotted. 
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